Critical Phase of Bond Percolations on Growing Networks 
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The critical phase of bond percolation on the random growing tree is examined. It is shown that 
the root cluster grows with the system size N as and the mean number of clusters with size s 
per node follows a power function n s oc s~ T in the whole range of open bond probability p. The 
exponent r and the fractal exponent ip are also derived as a function of p and the degree exponent 
7, and are found to satisfy the scaling relation r = 1 + ip -1 . Numerical results with several network 
sizes are quite well fitted by a finite size scaling for a wide range of p and 7, which gives a clear 
evidence for the existence of a critical phase. 
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I. INTRODUCTION 

The study of complex networks has been one of the 
most popular topics for many research fields in the last 
decade This activity has been drown by the dis- 

coveries of the small-world and the scale free (SF) @ 
properties, which are common to many real networks, 
e.g., Internet, WWW, social networks, food- webs. The 
former means the (sub-)logarithmic dependence of the 
mean shortest distance I between nodes with network 
size N as I oc log A (oc log NJ log log N) @, H H|, and 
the latter a power-law tail in the degree distribution 
P{k) oc fc -7 , where the degree k is the number of edges 
connected to a node. The dynamics on complex net- 
works, such as percolation, epidemic processes, interact- 
ing spin systems, coupled oscillators, have been exten- 
sively studied with stimulating our interests for the re- 
lationships between network topology and critical phe- 
nomena [7[. In most analytical approaches locally tree- 
like approximation is used to give many detailed physical 
pictures about critical phenomena on some network mod- 
els, particularly, on uncorrelated SF networks Q. On the 
other hand, many works have reported that systems on 
some gro wing networks show quite different phase tran- 
sitions d-[TJ. As shown in jTol . fllj , the percolation on 
some growing network models undergoes an infinite or- 
der transition with a BKT-like singularity; (i) the relative 
size of the giant component vanishes in an essentially sin- 
gular way at the transition point, so that the transition is 
of infinite order, and (ii) the mean number n s of clusters 
with size s per node (or the cluster size distribution in 
short) decays in a power-law fashion with s in the whole 
region where no giant component exists. Similar tran- 
sitions are observed for the interacting spin systems on 



a hierarchical SF network [12j | . and an inhomogeneous 
growing network [l4[ • 

The above unusual disordered phase is thought to be 
the same as the critical phase [TJI], which is also called 
patchy phase [l6[, observed on hierarchical SF networks 
and the Hanoi networks. By renormalization group (RG) 
techniques, Berker et al. [15| have studied bond perco- 
lations on the decorated (2,2)-flower [13], which is one 
of hierarchical SF networks, to show the existence of a 
critical phase, where RG flow converges onto the line of 
nontrivial stable fixed points. 

Let us turn our eyes upon exotic but regular graphs 
for a while. Critical phenomena on nonamenable graphs 
(NAGs), which arc defined as graphs with positive 
Cheeger constant, have been studied in recent years 
[l8l lisj . Roughly speaking, NAGs are regular graphs 
having small- world property I oc log N. Hyperbolic lat- 
tices and regular trees are typical examples of NAGs. 
It has been predicted that the bond percolation on a 
NAG exhibits a multiple phase transition (MPT) which 
takes three distinct phases according to the open bond 
probability p as follows; (i) the non-percolating phase 
(0 < p < Pel) in which only finite size clusters exist, 
(ii) the critical phase (p c i < p < p C 2) in which there 
are infinitely many infinite clusters, and (iii) the per- 
colating phase (p C 2 < p < 1) in which the system has a 
unique infinite cluster. Here infinite cluster means a clus- 
ter whose mass diverges with system size N as with 
< <p < 1. Recent paper [13, El (see also [HHU) per- 
formed Monte-Carlo simulations for the bond percolation 
on the enhanced binary tree, which is one of NAGs, to 
show that the system undergoes the MPT, and the crit- 
ical phase has the following properties; (i) continuously 
increasing of the fractal exponent ip with p, where ip is 
defined as 



so (AO oc N+, 
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so(N) being the mean size of cluster to which the root 
node belongs (hereafter we refer to this cluster as the root 



cluster), and (ii) the cluster size distribution n s 
having a power-law tail 

n s oc s~ r , 
with p-depcndent r satisfying 

r = 1 + ip' 1 . 



always 



(2) 



(3) 



We already know that standard systems on the Euclidean 
lattices (which are amenable) have just one critical point 
and no critical phase [24j . while it is unclear and rarely 
discussed so far whether percolations or other processes 
on complex networks have such a critical phase or not. 

The aim of this paper is to connect two concepts, un- 
usual phase transitions on growing networks and the 
MPT on NAGs. For this purpose we analyze the 
bond percolation on the growing random tree (GR tree) 
[25l [26j . As already mentioned, growing networks are 
expected to show a criticality in the region where its or- 
der parameter takes zero, while it is also trivial that any 
tree has no ordered phase due to the absence of loops. 
Thus, systems on the GR tree are always critical in the 
whole range of p (except p = 0,1). We show analytically 
that the exponents r and ip change continuously with the 
open bond probability p, with satisfying the scaling re- 
lation ((3]). We also perform the Monte-Carlo simulation 
to show that finite size scaling for n s is quite well fitted 
for any p, irrespective of the degree exponent 7. 



II. MODEL 

A realization Tjv of the GR tree with N nodes is ob- 
tained as a tree at time t = N generated by a stochastic 
process as follows: One starts at time t = 1 with Ti 
consisting of just one isolated node which we call the 
initial node or the root (for later convenience it is con- 
nected to a dangling bond so that the initial degree is 
k = 1). At each time step, a new node is added and 
linked to T^v to make Tjv + i (see figure [TJ. The probabil- 
ity that the new node is linked to a node of Tjy with the 
degree k is proportional to the linear attachment kernel 
Ak = k + a (a > —1). The stationary degree distribution 
P{k) of the resulting tree is known to have a power-law 
tail P(k) oc fc -7 , where the degree exponent 7 is related 
to a as 7 = 3 + a (2(| . 

We now consider the bond percolation on the GR tree 
with N nodes. Each bond is open with probability p, 
otherwise closed. 



III. DERIVATION OF ip AND r 

First, we evaluate the fractal exponent ip of the root 
cluster. Our calculation here is completely in the same 
way as for the Ising spin system [27[ ■ Let us assign a level 
to each node according to the distance I from the initial 
node. The initial node is regarded as at the 0-th level. 




FIG. 1: Example of growth of the growing random tree. 



Let nfj denote the mean number of nodes at the Z-th 



level on Tjy. Then the sum of all degrees of the nodes at 
the Z-th level is equal to + . ^ n adding a new 

node to Tjv, the probability that the new node is linked 



to any node at the Z-th level is ■ 



,(m) 

N 



)/[(2- 



a)N — 1], If this happens the new node itself is then at 
the (I + l)-th level, so we obtain 



,0+1) 
l N+l 



N 



(I) 



N 



C 2 N-1 



(I > 0), (4) 



where c\ = 1 + a, c 2 = 2 + a, and the boundary condition 



is riff = 1 for all N. The generating function of 11 



1=0 



(i) 

N ■ 



(5) 



is nothing but the mean size of the root cluster, and the 
recursion relation is obtained from (j4|) as 

(c 2 N - l)G N+1 (p) - (c 2 N + c lP )G N {p) - 1. (6) 

It is easily solved as 

1 



G N ( P ) 



1 + cip 

eg r(i-c 2 - 1 )r(^ + C2 - 1 Cl p) 

1 + cip r(l + c 2 - 1 c 1 p)T{N - C2" 1 ) ' 



(7) 



For N 3> 1, the second term becomes dominant, so that 
the mean size of the root cluster grows asymptotically as 



G N (p) ~ JV" 



l+(l + c)p 
N 2 + ° . 



and thus we find the fractal exponent ip as 
1 + (1 + a)p 1 + (7 - 2)35 



7-1 



(8) 



(9) 



Next, we derive the cluster size distribution n s . We 
consider the time evolution of the number N l s (t) of clus- 
ters which have s — 1 open bonds (s connecting nodes) 
and I closed bonds at time t. The time evolution of N l s (t) 
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FIG. 2: (Color online) Fractal exponent of the largest 
clusters (blue square) and ip of the root cluster (red circle) 
on the GR tree with a — 0,1, 100, from top to bottom. The 
dotted lines show the analytical prediction ((9} • iph and tp are 
given by the fit of those with several sizes. 



is given as 

2Vj(t + l) = N l s (t)- a ' S + l - 2 Nl(t) 



+P 

+q 



a'N - 1 

0/(8-1) + 1-2 , 



a'N - 1 



a ' S ,i l * N l -\t)+q5 s>1 5 lA , (10) 
a N — 1 



with the initial condition N s (t = 1) = 5 St ±6i t i, where 
a' = a + 2 and q = 1 — p. Assuming the asymptotic form 
as N l s (t) = tn l s = Nn l s (t > 1), we obtain 



(a's + I + a' - 2)n\ = p(a's + I - a' - 2)n l s _ 1 



(11) 



+q(a's + I - 3)n' s 1 + qa'6 sA 6i 
Summing over I gives the cluster size distribution as n s - 

EX; 



[a'ps + a' — 2p)n s = p(a s — a' — 2)n s -i 

I I 
which is evaluated if a — > oo as 



T[s + (l+p)/p] 



Hi OC S 



(12) 



(13) 



as previously reported in 0, [28| . To obtain an estimate 
of n s for finite offset a, we assume I* = E;^ n s/ n s for 
each s to satis fy s — 1 : I* = p : q. This approximation is 
the same as of [29|. Then n s is obtained as 

T[ S -(a'p + 2)/(a'p + q) + l] 
3 T[s + (a' -p-l)/(a'p + q) + lY { ' 

Thus its asymptotic form for s 3> 1 is given by @ with 

3 + p + ap + a 7 + (7- 2)p 



1 + p + up 



l + (7-2)p' 



(15) 
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FIG. 3: (Color online) Finite size scaling for the clus- 
ter size distribution n s of the GR tree with a — 0, and 
p = 0.3, 0.5, 0.7, from top to bottom. The number of nodes is 
taken 65536(red circle), 32768(blue square), 16384(green di- 
amond), 8192(black triangle), 4096(yellow inverted triangle). 
The dotted line follows the exponent of Eq. (|15[) . We obtain ip 
by numerical simulations, and determine r by scaling relation 
Q with numerically obtained 4>- 



Note that the scaling relation ^ is satisfied between ijj 
in Eq.© and r in Eq. (fT5")) . This result means that n s on 
the GR tree always follows power-law type, irrespective 
of 7 or p (0 < p < 1). 

We may say the same in terms of the correlation vol- 
ume and the correlation length. The connectedness func- 
tion d(l,p) of node i is defined as the probability that 
a randomly-chosen node, whose distance from node i is 
I, belongs to the same cluster with node i. Assuming 
that Ci(l,p) for p < p c i decays as a single exponential 
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function 

Ci{l,p)<xe- 1 ^, (16) 

where is the correlation length. The correlation volume 
of node i is given by summing the connected function 

as 

& = 52m(!)Ci(!, P ) t (17) 
i 

where rii(l) is the number of nodes whose distance from 
node i is I. t;Y is the size of cluster to which node i be- 
longs. Since rii(l) grows exponentially with I for graphs 
with the small-world property, it is possible that some 
correlation volumes diverge even when correlation length 
decays exponentially. As pointed in [20| , the critical 
points p c i and p C 2 of NAGs correspond to the points 
above which the correlation volume and length diverge, 
respectively. Similarly, these grow with p in a differ- 
ent way from each other on the GR trees: Eq.© in- 
dicates that the correlation volume of the initial node, 
$ , diverges if p > (p c i = 0), while the correlation 
length Q between any node pairs behaves like that of 
one-dimensional percolations, ~ — 1/lnp, and never 
diverges at any p < 1 (p C 2 = 1). 

Note that in the MPT there is a certain open probabil- 
ity p s between p c \ and p C 2 , above which the mean clus- 
ter size s oc J2 B s 2 n s diverges [3(|. Clearly, Eq. ((TS|) with 
t = 3 gives the value p s = (7 — 3)/2(7 — 2), so that p s de- 
pends on the degree exponent 7 in contrast to p c \ and p C 2- 
For 7 < 3, p s reaches in the limit N —> 00, so the mean 
cluster size diverges irrespective of p(> 0). Interestingly, 
another expression obtained by rewriting it in terms of 
the moments of k, p s = (7 - 3)/2( 7 - 2) = {k)/{k 2 - k), 
is also obtained for the percolation threshold on uncor- 
rected SF networks. A similar expression for the Ising 
case on the GR tree has already been obtained in [27| ■ 

IV. NUMERICAL RESULTS 

In this section, we check the above analytical results 
by numerical calculations. We perform the Monte-Carlo 
simulations for the bond percolation on the GR trees. 
The number of nodes is taken from about 4096 to 65536. 
To measure the mean size of the root clusters, we gen- 
erate 1000 graph samples from the GR tree and simu- 
late the bond percolation 5000 times on each sample. 
Figure [5] plots the fractal exponent of the root clus- 
ter ip and the analytical prediction Here we eval- 
uate ifj of the system with A nodes by the difference 
(In s (2 A) - In s (A/2))/(ln 2A-ln(A/2)). Numerically- 
obtained ip does not depend on A except for small values 
°f P {p ^ 0.2), where ip tends to approach the analytical 
line ([9]) with increasing A (not shown) . We also measure 
another critical exponent that characterizes the mean 
size of largest clusters s max (A) cx A'^ L . The estimated 
value in Fig[2] shows a good agreement with ip, which 



n s N* T 




FIG. 4: (Color online) Finite size scaling for the cluster 
size distribution n a of the GR tree with 01 — 1 (top panel), 
100 (bottom panel). Each curve takes open bond probability 
p = 0.2, 0.5, 0.8, from right to left. The number of nodes is 
taken 65536(red circle), 32768(blue square), 16384(green di- 
amond), 8192(black triangle), 4096(yellow inverted triangle). 
We obtain ijj by numerical simulations, and determine r by 
scaling relation ([3]) with numerically obtained if). 

indicates that the root clusters can be regarded as the 
largest ones, so(N) ~ s m ax(A)- To check n s , we assume 
a finite size scaling form [20| ; 

n s (N) = N-^fisN-f), (18) 

where the scaling function behaves as 

fl J rapidly decaying func. for 1 > 1, , . 
J\ x ) ~ I x -t f or x <^ 1 \ ) 

Here we also assume the scaling relation ^ . Figure [3] 
shows our finite size scaling for n s on the GR tree with 
7 = 3. Our scaling with several sizes are quite well fit- 
ted in a wide range of p. We also find similar results for 
various values of 7 as shown in FigJU We consider the 
scaling relation ([3]) is quite general. If n s oc s~ T holds 
asymptotically, a natural cutoff s max (A) of the cluster 
size distribution (a natural cutoff of the degree distribu- 
tion was introduced in [3ll ] ) is given as 

poo 

N / n s ds ~ 1 -> s max (A) oc A^r . (20) 
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Then Eq.Q follows with replacing ipi, by ip. 

Finally we note that a standard finite size scaling anal- 
ysis does not work to determine jp r i and p C 2 by using the 
data in the critical phase (see [2l| for the similar argu- 
ment of the enhanced binary tree). Suppose that in the 
critical phase s ma x oc holds with an increasing ana- 
lytic function of iph(p) for p c \ < p < p C 2- Then one can 
expand iph(p) around any p* in the critical phase so that 
■0L (p) — iph (j>* ) oc p — p* , which leads us to 



-V>l(p*) 



g(logN(p-p*)), 



(21) 



where g(.) is a "scaling function" around p* . This means 
that the scaled parameter p* is arbitrary as long as it 
is in the critical phase and its boundary p c \ and p c i is 
hardly determined from this type of analysis for practical 
use. 



critical behavior is similar to that observed in a NAG 
[20j . The critical phase is characterized by a power- law 
behavior of n s with varying exponent. The same prop- 
erty is observed on the decorated (2,2)-flower |l5[ and the 
Hanoi networks After the submission of this paper, 
Sato and the authors calculated ip and r of the deco- 
rated (2, 2) -flower by a generating functional approach 
to confirm that the phase is indeed critical everywhere 
[32j . As already mentioned, the critical phase is caused 
by opening a gap between the points at which the cor- 
relation volume and correlation length start to diverge. 
What is the structural factor of networks for such a sep- 
aration to occur? Finding the answer is, we believe, an 
essential step in understanding usual and unusual phase 
transitions on complex networks. 



V. SUMMARY 

Our results indicate that the bond percolation on the 
GR tree has a critical phase for < p < 1, and the 
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